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ABSTRACT 

We have made extensive spectroscopic and photometric observations of PG 1224+309, 
a close binary containing a DA white dwarf primary and an M4+ secondary. The Ha 
line is in emission due to irradiation of the M-star by the hot white dwarf and is seen 
to vary around the orbit. From the radial velocities of the Ha line we derive a period 
of P = 0.258689 ± 0.000004 days and a semi-amplitude of K Ha = 160 ± 8 km s _1 . We 
estimate a correction AK = 21 ± 2 km s _1 , where Km = K# a + AK. Radial velocity 
variations of the white dwarf reveal a semi-amplitude of Awd = 112 ± 14 km s _1 . The 
blue spectrum of the white dwarf is well fit by a synthetic spectrum having T e g = 29, 300 
K and logg = 7.38. The white dwarf contributes 97% of the light at 4500 A and virtu- 
ally all of the light blueward of 3800 A. No eclipses are observed. The mass inferred for 
the white dwarf depends on the assumed mass of the thin residual hydrogen envelope: 
0.40 < M WD < 0.45 M for hydrogen envelope masses of < M H < 4 x 10" 4 M . We 
argue that the mass of the white dwarf is closer to 0.45 M , hence it appears that the 
white dwarf has a relatively large residual hydrogen envelope. The mass of the M-star is 
then M M = 0.28±0.05 M , and the inclination is i = 77±7°. We discuss briefly how PG 
1224+309 may be used to constrain theories of close binary star evolution, and the past 
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and future histories of PG 1224+309 itself. The star is both a "post-common envelope" 
star and a "pre-cataclysmic binary" star. Mass transfer by Roche-lobe overflow should 
commence in about 10 10 yr. 

Subject headings: stars: binaries: close — stars: individual (PG 1224+309) — stars: 
variables — stars: white dwarfs 
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1. Introduction 

PG 1224+309 was cataloged as a UV-excess object 
in the Palomar-Green (PG) survey (Green, Schmidt, 
&: Liebert 1986). It also appears in the Tonantz- 
intla survey as Ton 617 (Iriarte & Chavira 1957; 
Chavira 1958) and the Case low-dispersion northern 
sky survey as CBS 60 (Sanduleak & Pesch 1984). 
Its mean magnitude and colors are V = 16.164, 
B — V = -0.065, V - R = 0.038, and R - I = 0.312 
(this work). This object was studied by Ferguson, 
Green, & Liebert (1984) who were searching for cat- 
aclysmic variable-like stars in the PG survey. Fer- 
guson et al. (1984) showed that PG 1224+309 was a 
binary, identified the hot star as a DA white dwarf 
with T eff = 28, 000 K and log g = 7.7 ± 0.8, and clas- 
sified the cool star as an M2 dwarf based on its col- 
ors. Refined values of T cff = 29, 300 ± 1000 K and 
\ogg = 7.38 ± 0.1 for the white dwarf were later 
provided by J. W. Liebert & P. Bergeron (1997). 
PG1224+309 was not discussed again in the liter- 
ature for thirteen years. Orosz, Wade, & Harlow 
(1997) looked for radial velocity variations in a num- 
ber of stars that were characterized as composite spec- 
trum binaries by Ferguson et al. (1984), and found 
for PG 1224+309 that the radial velocity of the M 
star changed by 83 + 11 km s^ 1 between two obser- 
vations w 0.9 days apart. Orosz et al. (1997) also 
showed that PG 1224+309 has a strong and variable 
Ha emission line. Photometric observations by S. Bell 
(private communication 1996) showed that the light 
varies smoothly on a timescale of hours. 

All of these observations strongly suggested that 
PG 1224+309 is a close binary system, in which the M 
star is irradiated by the white dwarf. We therefore ob- 
tained additional photometric and spectroscopic data 
on this object, reported in this paper. Our goals were 
to characterize the orbit, to estimate the masses and 
other properties of the two stars, and to place the sys- 
tem in its evolutionary context. The additional data 
confirm that PG 1224+309 is a short-period, double- 
lined spectroscopic binary system with photometric 
and emission-line variations evidently caused by illu- 
mination of one hemisphere of the cool M star by the 
hot white dwarf. We present below a description of 
our observations and data analysis. We discuss vari- 
ous constraints on the component masses and the ge- 
ometry of the system and end with a brief discussion 
of the evolutionary status of this binary star. Systems 
such as PG 1224+309 may illuminate difficult aspects 



of common-envelope binary evolution. 
2. Observations 

2.1. MDM^j Spectroscopic Observations 

We obtained 55 exposures between 1996 Decem- 
ber and 1998 March using either the 2.4 m Hiltner 
Telescope and modular spectrograph, or the 1.3 m 
McGraw-Hill Telescope and Mark III spectrograph. 
Table | summarizes the observations. The 2.4 m data 
had a spectral resolution of ~ 3.5 A FWHM, and the 
1.3 m data had ~ 5 A, as determined from compar- 
ison lines. All spectra covered H/3 to Ha, and the 
later 2.4 m spectra (taken with the 2048 x 2048 CCD) 
covered from 4000 to 7500 A with considerable vi- 
gnetting toward the ends. Individual exposures were 
300 to 600 s at the 2.4 m telescope, and typically 900 
s at the 1.3 m. All stellar spectra were bracketed 
in time with comparison-lamp exposures, and typical 
rms residuals of the wavelength solutions were < 0.1 
A. Standard IRAF^] specred tasks were used to sub- 
tract the electronic bias, divide by fiat-fields, opti- 
mally extract one-dimensional spectra, subtract sky 
background, and rebin the spectra onto a uniform 
wavelength scale. We obtained spectra of hot stars 
and flux standards to convert our spectra to absolute 
flux units. 

At the 2.4 m we used a narrow slit (1 arcsec), so our 
spectrophotometry suffers from variable, uncalibrated 
losses at the slit jaws; we estimate that this effect in- 
troduces a ~ 30 percent uncertainty in our absolute 
flux levels. In addition, continuum shapes derived 
with the modular spectrograph sometimes suffer 10- 
20 per cent errors in their overall slopes, which we do 
not understand. We did rotate the spectrograph to 
align the slit with the parallactic angle whenever the 
zenith distance was substantial, so differential refrac- 
tion is not the culprit. 

2.2. Kitt Peak Photometric Observations 

We obtained images of PG 1224+309 during the in- 
terval 1998 May 5-13 with the KPNO 0.9m telescope, 
which was equipped with the T2KA 2048 x 2048 CCD 
at both the f/ 13.5 focus (May 5-8) and the f/7 focus 



5 MDM Observatory is operated by a consortium of Dartmouth 
College, Columbia University, the Ohio State University, and 
the University of Michigan. 

6 IRAF is distributed by the National Optical Astronomy 
Observatories. 



3 



Table 1 

Log of Spectroscopic Observations 



Start a Finish a N Setup Detector 



447.03 


453.05 


8 


MDM 2.4+mod 


Tekl024 


486.89 


491.96 


18 


MDM 1.3+MkIII 


Tekl024 


494.94 


498.90 


8 


MDM 2.4+mod 


Tekl024 


623.66 


630.66 


7 


MDM 2.4+mod 


SiTe2048 


800.06 


801.06 


2 


MDM 2.4+mod 


SiTe2048 


842.04 


846.98 


8 


MDM 2.4+mod 


SiTe2048 


894.90 


895.98 


4 


MDM 2.4+mod 


SiTe2048 


955.83 


955.91 


8 


Lick 3m + Kast 


Reticon 1200 x 400 


987.76 


987.84 


8 


Lick 3m + Kast 


Reticon 1200 x 400 



a Heliocentric Julian date minus 2,450,000. 



(May 9-13). Typical exposure times were 400 s for 
the B filter, 300 s for V, 200 s for (Kron-Cousins) 
R, and 250 s for (Kron-Cousins) /. The night of May 
4/5 was photometric and we obtained 24 observations 
of 16 different standard stars from the list of Landolt 
(1992) in each of the B, V, R, and / filters, in ad- 
dition to the nine B, V, R, and I sequences on PG 
1224+309. 

Standard IRAF tasks were used to subtract the 
electronic bias and perform flat-field corrections. We 
used the programs DAOPHOT He, ALLSTAR and 
DAOMASTER (Stetson 1987; Stetson, Davis, & Crab- 
tree 1991; Stetson 1992a, b) to compute the photomet- 
ric time series of PG 1224+309 and five nearby field 
stars. The observations of the Landolt standard stars 
were used to derive the transformation from DAO- 
MASTER instrumental magnitudes to the standard 
scales. The formal errors for the zeropoints of the 
calibrated magnitude scales are 0.018 mag for the B 
filter, 0.005 mag for the V filter, 0.012 mag for the 
R filter, and 0.013 mag for the I filter. Another in- 
dication of the quality of the photometric calibration 
comes from an examination of the measured compar- 
ison star magnitudes. Each comparison star had nine 
measures in the four different filters. There was very 
little scatter in these measurements: the standard de- 
viations ranged from 0.003 to 0.030 mag. 



2.3. Lick Spectroscopic Observations 

Blue spectra of PG 1224 + 309 were obtained 
with the 3m Shane Telescope and Kast spectrograph 
(Miller & Stone 1993) at Lick Observatory on 1998 
May 22 UT and again on 1998 June 23 UT. On each 
occasion the star was observed for approximately 2 
hours during which a set of eight 15-minute exposures 
was obtained, covering about 1/3 of the orbital cycle. 
The purpose was to measure the radial velocity curve 
of the white dwarf from the high-order Balmer lines 
in its spectrum, so the observations were planned to 
sample a different part of the orbital cycle of the sys- 
tem on each run. The spectra cover the range 3200- 
4550 A with a spectral resolution of 2.8 A. They were 
obtained through a 2" slit using an 830 mm -1 grism. 
The seeing varied between 1"5 and 2" over the course 
of the observations resulting in a variable loss of light 
at the slit jaws. Moreover, since the star had to be 
followed through a fairly high airmass, some spectra 
suffer from additional light losses at the shortest wave- 
lengths because of differential atmospheric refraction, 
in spite of our efforts to keep the slit at the parallactic 
angle. A log of the observations is included in Table 

& 

The spectra were reduced in a standard manner, as 
described earlier, although we used our own reduction 
software while making only limited use of the IRAF 
software package. In the final individual spectra, the 
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Table 2 
Ha Radial Velocities 



HJD a 


V 


HJD a 


V 


HJD a 


V 


HJD a 


V 




(km s" 1 ) 




(km s- 1 ) 




(km s _1 ) 




(km s" 1 ) 


447.998 


-48 


490.777 


-73 


624.656 


57 


846.017 


163 


451.052 


125 


491.953 


2 


624.660 


47 


846.025 


161 


451.058 


120 


491.964 


-32 


626.658 


130 


846.962 


-113 


487.032 


20 


497.888 


86 


626.662 


172 


846.969 


-129 


487.833 


-112 


497.894 


77 


630.658 


-156 


895.968 


144 


488.859 





498.894 


187 


800.062 


14 


895.976 


138 


488.870 


9 


498.900 


148 


842.038 


-154 






489.065 


209 


623.658 


-75 


845.010 


173 






489.840 


215 


623.665 


-97 


846.010 


147 







a Heliocentric JD of mid- integration minus 2450000. 



signal-to-noise ratio at 4200 A ranges between 15 and 
25. The wavelength scale was originally derived from 
spectra of arc lamps obtained at the beginning of each 
sequence of object spectra. The wavelength scales 
of individual spectra were then refined by using the 
emission-line spectra of the night sky recorded during 
the same exposure to measure small relative offsets 
and rectifying them. As a result of this procedure, 
the relative velocity scales of individual spectra agree 
to 8 km s _1 or better. 

3. Data Analysis 

3.1. The Spectroscopic Period and the Ha 
Emission Line Radial Velocity Curve 

We used the set of spectra from MDM to deter- 
mine the spectroscopic period since the time cover- 
age for these data is by far the longest. The MDM 
spectra (discussed later) showed a narrow, variable 
Ha emission line superposed on a strong white dwarf 
spectrum of type DA (see also Orosz, Wade, & Har- 
low 1997). For 33 spectra with adequate signal-to- 
noise, we fitted this feature with a Gaussian, adjust- 
ing approximately for the sloping baseline caused by 
the white dwarf absorption. Table || gives the result- 
ing radial velocities. A Gaussian-convolution algo- 
rithm (Schneider & Young 1980) gave similar results, 
but with slightly greater scatter around the best fit 
(below). We searched for periods by constructing a 



highly oversampled grid of trial frequencies spanning 
to 8 cycles d _1 , fitting a least-squares sinusoid to 
the velocities at each frequency, and plotting the in- 
verse of the fit variance as a function of frequency 
to form a 'residual-gram' (Thorstensen et al. 1996). 
This method works well when the modulation is si- 
nusoidal and its amplitude exceeds the observational 
uncertainty. Figure |l| shows the result, which is con- 
sistent with a single frequency near 3.86 d _1 ; the cycle 
count is unambiguous. A least-squares sinusoid fit of 
the form 



v(t) = 7H Q + Kji a sin 



2ir(t - To) 



P 



gave 



T 
P 

7h q 
a 



HJD 2450624.537 ±0.002, 
0.258689 ± 0.000004 d, 
165 ±9 km s" 1 , 
10 ± 6 km s _1 , and 
29 km s" 1 , 



where a is the uncertainty of a single measurement 
inferred from the goodness-of-fit, and the uncertain- 
ties are formal 1-tr errors. Discarding the two worst- 
fitting points changed K, 7, and a to 160 ±8, 8 ± 5, 
and 26 km s _1 respectively, but did not significantly 
affect the ephemeris. Figure § (lower panel) shows all 
the velocities folded on the ephemeris. The sinusoid 
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PG1 224+309 Period Search 
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JuLjulitJiiJ 



10 
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Fig. 1. — Results of a period search on the Ha emis- 
sion velocities. The vertical axis is the inverse of 
the variance around a sinusoidal fit at each trial fre- 
quency. In order to compress the large number of 
data points we show local maxima in this function, 
joined by straight lines. 

shown has been fitted to all but the two worst points. 
We adopt K Ka = 160 ± 8 km s" 1 . 

The Ha emission fits also produced estimates of 
the equivalent width (EW). These were noisier than 
the radial velocities, but showed a strong modulation 
at the same period. A sinusoidal fit to the EW time 
series, with the period fixed at the more accurately- 
determined radial velocity period, showed that the 
EW modulation lags in phase behind the radial ve- 
locities by 0.256 ± 0.016 cycle. This is nicely con- 
sistent with the 0.25-cycle lag expected if the emis- 
sion line arises on the side of the normal star fac- 
ing the white dwarf (Thorstensen et al. 1978; Venncs 
& Thorstensen 1994). The upper panel of Figure ^ 
shows the equivalent widths folded on the best period, 
together with a sinusoidal fit with a half-amplitude of 
3.0 ± 0.3 A and a mean of 4.4 ±0.2 A (formal er- 
ror). The sinusoidal fit is only fiducial, as the mod- 
ulation may not be strictly sinusoidal (Thorstensen 
et al. 1978), and our EW accuracy and sampling are 
not sufficient for us to detect small deviations from 
a sinusoid. In any case, the deep modulation of the 
equivalent widths, and the accurate agreement of the 
phase offset with the illumination model, make it very 
likely that the bulk of the Ha emission arises from flu- 
orescence on the illuminated face of the normal star. 
It is possible that the M dwarf also has intrinsic Ha 
emission, but it apparently does not dominate. 
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Fig. 2. — Lower panel: radial velocities of the Ha 
emission lines, derived from Gaussian fits to the line 
profile, folded on the best period, with the best-fitting 
sinusoid superposed. The points represented by open 
squares were excluded from the fit. Upper panel: 
Equivalent widths of the Ha emission line folded on 
the best period. The fiducial sinusoid shown is shifted 
by 0.25 cycle from the velocity ephemeris. 

3.2. Photometric Variations 

The B, V, R, and / light curves of PG 1224+304 
display approximately sinusoidal modulations with a 
period near « 6 hours, as expected based on the spec- 
troscopic analysis. We fit four-parameter sinusoids 
to each of the four light curves (using the spectro- 
scopic period P sp0 ct as the initial period guess) to 
determine the photometric period, semi-amplitude, 
mean magnitude, and the time of photometric min- 
imum. The results are displayed in Table ^[ The 
errors on the individual magnitudes were scaled to 
give xt = 1 f° r each curve and the errors on the fit- 
ted parameters were computed from the scaled er- 
rors. The errors listed in Table || may have been 
slightly underestimated, since the light curve shapes 
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Table 3 
Photometric parameters 



Filter 


N 


Period 


Semi-amplitude 


Mean 


T (photo) 


Spectroscopic phase 






(days) 


(mag) 




(HJD 2450000±) 


of T (photo) 


B 


38 


0.2587 ±0.0003 


0.034 ± 0.002 


16.099 ±0.002 


938.840 ±0.004 


1214.98 ±0.03 


V 


55 


0.2588 ±0.0001 


0.070 ±0.002 


16.164 ±0.001 


938.845 ± 0.002 


1215.00 ±0.02 


R 


39 


0.2586 ±0.0001 


0.104 ±0.003 


16.126 ±0.002 


938.847 ±0.002 


1215.01 ±0.02 


I 


50 


0.2582 ±0.0001 


0.098 ±0.002 


15.812 ±0.001 


938.849 ±0.001 


1215.02 ±0.02 



j r ^ 





Fig. 3.— The B, V, R, and J light curves folded 
on the adopted (mean) cphcmcris arc shown in the 
left panels, where each point has been plotted twice. 
The right panels show the folded light curves binned 
into ten phase bins, where the error bars represent 
the errors of the mean in each bin. Sine curves have 
been added to demonstrate the absence of eclipses; 
the minimum of each curve occurs at phase zero, and 
the semiamplitudes in the B, V, R, and I bands are 
0.04, 0.07, 0.10, and 0.10 mag respectively. 



may not be exact sinusoids. The periods derived for 
the four light curves are close to the spectroscopic 
period. We also computed the spectroscopic phase 
of the time of the photometric minimum for each 
light curve. All of these phases are consistent with 
a whole number of cycles, as expected if the cause 
of the modulation is due to irradiation of the M star 
by the white dwarf. We used the weighted mean of 
the photometric To and the spectroscopic To (prop- 
agated forward by 1207 cycles) to define the phase 
zero-point for the adopted ephemeris: To (mean) = 
HJD 2, 450, 938.8430 ± 0.0026. Our adopted orbital 
parameters are given in Table H, and the folded light 
curves are displayed in Figure || 

3.3. The White Dwarf Radial Velocity Curve 

The radial velocity of the white dwarf was mea- 
sured from each Lick spectrum by comparing it to a 
synthetic white dwarf spectrum, constructed with the 
parameters reported by Liebert & Bergeron (1997). 
The method involves assuming a radial velocity for 
the white dwarf, applying a Doppler shift and appro- 
priate scaling to the model spectrum, and comparing 
it to the observed spectrum in the range 3890-4400 A 
by means of the \ 2 test. By scanning a range of pos- 
sible white dwarf radial velocities in very fine steps 
we obtain the best estimate as the one with the low- 
est value of x 2 as well as error bars corresponding to 
68% (Tct) confidence limits. The measured heliocen- 
tric radial velocities are given in Table |5| along with 
their corresponding orbital phases, computed from 
the adopted mean ephemeris given in Table In Fig- 
ure [| we plot the radial velocity of the white dwarf as 
a function of this orbital phase. 

The observed radial velocities were fitted with a 
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Table 4 
Adopted Orbital Parameters 



parameters value 

Pspcct (days) 0.258689 ± 0.000004 

K Ha (km s^ 1 ) 160 ±8 

A WD (km s" 1 ) 112 ±14 

To(spect) (HJD 2,450,000+) 626.608 ± 0.002 

T (photo) (HJD 2,450,000+) 938.844+0.003 

T (mean) (HJD 2,450,000+) 938.843 ± 0.003 



Note. — To refers to the sine-fit of radial velocities 
measured from the Ha emission line, i.e. inferior con- 
junction of the M star. 



PG 1224 + 309: W.D. radial velocity curve 




0.3 0.4 0.6 0.8 1 

Orbital Phase (cycles) 

Fig. 4. — The measured radial velocity curve of the 
white dwarf with two circular-orbit models super- 
posed for comparison. The radial velocities were mea- 
sured from the high-order Balmer lines (Hy, H5, and 
He) and have been corrected for the Earth's motion 
around the Sun. The orbital phase is computed ac- 
cording to the adopted (mean) ephemeris. The solid 
line shows a fit in which the phase of inferior conjunc- 
tion of the white dwarf was treated as a free parame- 
ter, while the dashed line shows a model in which this 
phase was fixed at 0.5. 

circular-orbit model. Because the observations do not 
cover a complete orbital cycle we cannot use them 
to obtain a measurement of the period. Hence we 
adopt the period derived from the Ha velocities and 
we express the model as 



v{4>) = 7wd + Kwd sin [27r(0 - Acf>)] 



(1) 



where <p is the orbital phase from the adopted ephemeris 
(0 < <f> < 1) and A0 is the phase of inferior conjunc- 
tion of the white dwarf, which is expected a priori to 
be 0.5. The free parameters of the model are 7wd> 
A'wd and Ac/). A minimum \ 2 fit of the model to the 
data gives the following parameters: 



A0 
A'wd 
7wd 



- 442+°' 049 
= 119±|j km s- 1 
= 5 ±21 kins" 1 



and 



with errors corresponding to 68% confidence limits. 
Since the best-fitting value of Acf> is consistent with 
expectation (according to the F-test), we may fix A</> 
to its expected value. Under this assumption we ob- 
tain: 



Ac/) 
7wd 



0.5 (fixed), 

112 ±14 km s-\ 

-4 ±12 km s- 1 



and 



The best-fitting model radial velocity curves (with A</> 
both free and fixed) are superposed on the data in 
Figure ^, for comparison. 

3.4. White Dwarf Atmospheric Parameters 

Using the fitted radial velocity curve of the white 
dwarf, we corrected the Doppler shifts in the Lick 
spectra and averaged them together to obtain a spec- 
trum with a high signal-to-noise ratio, which is free 
from orbital Doppler broadening. We chose to use 
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Table 5 

Heliocentric Radial Velocities from White Dwarf Absorption Lines 



HJD a 




v G 


HJD a 










(km s x ) 






(km s 1 ) 


955.82971 


0.665 


38 ±49 


987.75956 


0.094 


-152 ±40 


955.84151 


0.710 


116 ±41 


987.77064 


0.137 


-160 ±45 


955.85258 


0.753 


130 ± 42 


987.78170 


0.180 


-96 ± 44 


955.86365 


0.796 


104 ± 42 


987.79276 


0.222 


-65 ± 46 


955.87471 


0.839 


122 ± 46 


987.80385 


0.265 


-114 ±50 


955.88578 


0.881 


56 ±59 


987.81493 


0.308 


-47 ±57 


955.89683 


0.924 


30 ±56 


987.82601 


0.351 


-55 ±56 


955.90792 


0.967 


35 ±53 


987.83709 


0.394 


30 ±78 



a Heliocentric JD of mid- integration minus 2450000. 

b Orbital phase, computed from the adopted (mean) ephemcris. 



four of the Lick spectra for this purpose, namely the 
ones spanning orbital phases 0.09 to 0.22. These four 
spectra were obtained under the best seeing condi- 
tions and at the lowest airmass. As a result they in- 
dividually have a high signal-to-noise ratio and the ef- 
fects of differential atmospheric refraction are negligi- 
ble. Also, the line-emitting face of the secondary star 
is directed mainly away from the observer at these or- 
bital phases. Hence the contamination of the cores of 
the absorption lines in the white dwarf spectrum by 
narrow emission lines is minimal. The average spec- 
trum of the white dwarf is shown in Figure |5J. 

We used this Lick "restframe" spectrum with its 
improved signal-to-noise to estimate values of the 
white dwarf surface gravity \ogg and effective tem- 
perature T cS . (Recall that T cff = 29,300 K and 
\ogg — 7.38 ± 0.1 according to Liebert & Bergeron 
[1997]. Their uncertainties include estimates of the 
systematic errors due to modeling assumptions and 
input physics.) We computed a grid of pure hydrogen, 
LTE model spectra, covering the range log g = 7.18 
to 7.58 in steps of 0.1 and the range T c s = 28,300 
K to 29,800 K in steps of 500 K. These synthetic 
spectra were computed using I. Hubeny's programs 
TLUSTY and SYNSPEC, and the broadening of the 
hydrogen lines was treated using the tables of Schoen- 
ing and Butler (see Hubeny & Lanz 1997). We fitted 
the Lick restframe spectrum over the wavelength in- 



terval 3700-4533 A, excluding the Ca II H and K lines 
and features at approximately 4018 and 4174 A. 

The observed spectrum of PG 1224+309 consists 
of three components: the white dwarf spectrum, the 
spectrum of the late- type secondary (see below), and 
the spectrum from the heated face of the secondary. It 
is likely that the unilluminated portion of the M star 
does not make a noticeable contribution in the U and 
B bands, but the same thing cannot be said about the 
heated face of the M star. The Lick restframe spec- 
trum was constructed so that the contamination from 
the heated face was minimized as much as possible, 
but the fact that the Ha emission line is always visible 
suggests that some of the heated hemisphere is visible 
near phase 0.0. To allow for this extra light, we used 
a fitting method similar to the one given in Marsh, 
Robinson, & Wood (1994), which can be summarized 
as follows. The data and models are all normalized 
at a common wavelength (4501 A). Each model spec- 
trum is scaled by an amount w (0 < w < 1) and 
subtracted from the data, a third-order polynomial is 
fitted to the residuals, and the rms residual of the fit 
is recorded. The optimal fit is defined to be the one 
that gives the lowest overall value for the rms, cor- 
responding to the "smoothest" difference spectrum. 
Here we are assuming that the spectrum of the irra- 
diated hemisphere is featureless. The best decompo- 
sition is shown in Figure 0. The model white dwarf 
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spectrum has T e g = 29,300 K and logg = 7.38 and 
contributes 97% of the light at 4500 A and virtually 
all of the light blueward of 3800 A. Our values of T e g 
and log g are the same within the errors as those found 
by Liebert & Bergeron. Thus we adopt logg = 7.38 
for the discussion below. 




3800 4000 4200 4400 

Wavelength (A) 

Fig. 5. — Upper panel, from top to bottom: The blue 
"restframe" spectrum of PG 1224+309, constructed 
from four spectra having a mean orbital phase of 
4> — 0.16 (filled circles); the best-fitting model spec- 
trum (with \ogg = 7.38 and T cfi = 29,300 K), scaled 
by 0.967 (smooth line); and the difference spectrum 
with error bars on the individual points, offset up- 
wards by 0.4 for clarity. The horizontal bars with 
gaps indicate the wavelength regions used in the fit. 
Lower panel: The best-fitting model (smooth line) is 
shown with the detrended data, made by subtracting 
the polynomial fit to the residuals from the observed 
spectrum. 



3.5. The Spectral Type of the Secondary Star 

Having knowledge of the orbital ephcmeris, we 
searched in the red spectra for the contribution of the 
M star and for any variations of its spectrum with 
orbital phase. Figure |] shows our large-format MDM 
2.4 m spectra (which cover to A = 7500 A) averaged 
into 1/4-cycle phase bins centered on the four cardi- 
nal phases (conjunction and quadrature). Unfortu- 
nately, only a single spectrum fell in the bin nearest 
inferior conjunction of the M star. Nonetheless, a 
heating effect is clearly seen. In the lowermost spec- 
trum (inferior conjunction), the characteristic bumpy 
continuum of an M dwarf is strongly present; it is less 
pronounced, though quite obvious, in the quadrature 
phases. However, in the spectrum showing the largest 
portion of the illuminated face (third from bottom), 
the M-dwarf contribution is not obvious despite good 
signal-to-noise. Evidently the spectrum of the illumi- 
nated face is of earlier type than the unilluminated 
face (or is otherwise lacking strong molecular bands), 
so that the bands in the observed spectrum are greatly 
diluted. The illumination of the M star is discussed 
further below. 

To quantify the spectral contribution of the sec- 
ondary star and estimate its type, we prepared an av- 
erage of the eight large-format spectra taken within 
0.2 cycle of inferior conjunction. We then systemati- 
cally scaled and subtracted spectra from a library of 
M dwarfs. These stars were observed with the same 
instrumental setup, and have classifications from Boe- 
shaar (1976). We examined the subtracted spectra 
looking for good cancellation of the M-dwarf features. 
Types M3 and earlier tended to leave a strong over- 
subtraction of the Na I D lines when the scale factor 
was set to cancel the bands in the 6500-7500 A re- 
gion, and the bands themselves did not cancel well. 
A spectrum of type M5 did not subtract satisfacto- 
rily. The only acceptable results were from two stars 
classified as M4+ (Boeshaar 1976) or M4.5 (Henry, 
Kirkpatrick, & Simons 1994). Figure [?] shows the 
best decomposition. We conclude that the intrinsic 
spectrum of the M dwarf is most probably a little 
later than M4, and that the flux contributed at this 
phase by the unilluminated part of the M dwarf is 
F A ~ 1.5 x 10~ 16 erg cm" 2 s" 1 A" 1 at 6500 A, or 
15-20 percent of the light at that wavelength. The 
calibration uncertainties discussed earlier should not 
seriously affect these conclusions, since the features 
cover a fairly limited wavelength range. For compar- 
ison, Orosz, Wade, & Harlow (1997) adopted a spec- 
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5000 5500 6000 6500 7000 
Wavelength [A] 

Fig. 6.— Average spectra of PG 1224+309 in 0.25- 
cycle phase bins. Each is labeled by its mean phase. 
Only one spectrum fell in the bin nearest zero phase. 
No scale factors are applied, and the spectra are off- 
set from each other vertically by one unit (10~ 15 erg 
cm~ 2 s _1 A^ 1 ). 

tral type of M1-M3 with a contribution of 10% over 
the region 5320-6532 A, based on two low signal-to- 
noise spectra obtained near orbital phases 0.89 and 
0.53. 

Very few features are seen in the spectrum aside 
from the contributions discussed above. There is weak 
emission present at A5169, apparently also modulated 
by the illumination phase. This feature is fairly of- 
ten seen in cataclysmic binaries, where it is usually 
attributed to Fc II (Taylor & Thorstcnscn 1996). Cu- 
riously, there is very little He I emission; He I is 
conspicuous (though weaker than the Balmcr lines) 
in the illumination-effect stars Feige 24 (Thorstenscn 
et al. 1978) and EUVE 2013+400 (Thorstcnscn ct 
al. 1994). 

Because the Ha emission line arises mostly from 
one hemisphere of the M dwarf, its velocity ampli- 
tude is likely to underestimate the star's center-of- 




Wavelength [A] 



Fig. 7. — Decomposition of the spectrum. The top 
trace is the average of eight spectra taken within 0.2 
cycles of inferior conjunction of the M dwarf. The 
lower trace is a scaled spectrum of Gliese 896B, spec- 
tral type M4+ (Boeshaar 1976) or M4.5V (Henry, 
Kirkpatrick, & Simons 1994). The middle trace is 
the difference of the two. The top trace is offset up- 
ward by one unit (10~ 15 erg cm~ 2 s^ 1 A -1 ) to avoid 
overlap with the middle trace. 

mass motion. To account for this effect we tried to 
measure velocities of the M dwarf by cross-correlation 
methods, using our M dwarf library spectra (some of 
which have very accurate velocities from Marcy et 
al. 1987) as templates. Unfortunately, the signal-to- 
noise ratio of the available spectra proved inadequate 
for this purpose. 

4. Discussion 

4.1. Mass Constraints 

Knowing the surface gravity of the white dwarf in 
PG 1224+309 enables us to estimate its mass. The 
mass-radius relation for white dwarf stars does de- 
pend on the composition of the white dwarf and its 
position along the cooling track. In addition, the 
radius will be larger for a white dwarf that has a 
hydrogen layer atop the degenerate core. Liebert & 
Bergeron (1997) estimated the mass of PG 1224+309 
to be Mwd = 0.36 Mq based on their estimates 
of T e ff = 29,300 and logg = 7.38 and using then- 
available evolutionary tracks for carbon-oxygen white 
dwarfs (cf. Wood & Winget 1989). For these same 
values of log g and T e g, evolutionary tracks for he- 
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Hum white dwarfs by Althaus & Benvenuto (1997) 
yield M W d = 0.40 Mq (their Fig. 13). These mod- 
els do not include a hydrogen envelope. Benvenuto 
& Althaus (1998) published additional cooling tracks 
for helium white dwarfs that include a thin hydro- 
gen envelope. At log g = 7.38, these models im- 
ply 0.40 < Mwd < 0.45 Mq for hydrogen envelope 
masses of < Mh < 4 x 10~ 4 Mq. For the discussion 
below we will adopt 0.35 < M WD < 0.45 Mq. 
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Fig. 8. — Top: The mass of the M-dwarf secondary 
star as a function of AK for Mwd = 0.35 M , 
M WD = 0.40 Mq, and M WD = 0.45 Mq, assum- 
ing K WB = 112 ± 14 km s" 1 and K Ha = 160 ± 8 
km s _1 . The sizes of the statistical errors on Mm 
are illustrated by the three points with error bars. 
Bottom: The derived inclination as a function of 
AK for M W d = 0.35 M Q , and M W d = 0.40 M Q , 
Mwd = 0.45 Mq. The dashed lines indicate the range 
of allowed solutions for the Mwd = 0.40 Mq case, 
propagating one-sigma errors in ifwD and K^ a - 



The measured orbital velocities of the two stars de- 
termine the mass ratio: Q = Mm/Mwd = ^wd/^m- 
If Q were known accurately, then Mm and the to- 
tal mass of PG 1224+309 would simply scale from 
Mwd as estimated above. The velocity semiampli- 
tude of the white dwarf was measured reliably from 
the higher Balmer lines visible in the Lick spectra: 
Kwd = 112 ± 14 km s _1 . On the other hand, 
the M star's orbital speed is not directly measured, 
since the center of light of the Ha line emission docs 
not correspond to the center of mass of the M star, 
Km = Kft a + AK. Since the Ha emission line is 
formed in the heated face of the M star, we expect 
AK > 0, or Q = K wr >/(K Ha + AK) < 0.700 ±0.094, 
using the adopted Kn a = 160 ± 8 km s _1 . 

We can express AK as 

AK= — (1 + Q)K M 
a 

where AR is the shift of the "effective center" of the 
Ha emission line with respect to the center of mass of 
the M star (Wade & Home 1988) and the orbital sepa- 
ration is a = 1.71(Mtotai/M0) 1 / 3 Rq. In the case of a 
synchronously rotating M star with uniform Ha emis- 
sion over one hemisphere and no Ha emission over 
the other, AR = 4Rm/(3tt), where Rm is the radius 
of the M star. To complete an estimate of AK, we 
temporarily adopt values for R M of 0.25 ± 0.01 Rq 
and for M M of 0.25 ± 0.01 M , based on the dM4e 
eclipsing binary CM Draconis (Lacy 1977). Using 
Mwd = 0.4 ± 0.05 M Q , we then find Q = 0.625 and 
Mtotai = 0.65 Mq. Using K Ha = 160 ± 8 km s 1 , we 
find AK = 21 ± 2 km s" 1 . 

The true value of AK probably cannot be lower 
than our rough estimate of 21 km s , since AR will 
be larger if the Ha emission is more concentrated to 
the regions of the M star closest to the white dwarf 
(we argued above in Section 3.1 that the M star prob- 
ably has very little intrinsic Ha emission). In prin- 
ciple, AR could be as large as Rm {AK = 49 km 
s _1 ), but then one would not expect to see the equiv- 
alent width of Ha varying sinusoidally, as actually ob- 
served. Tighter upper limits on Q and AK come from 
considering the orbital inclination, which must be less 
than about 80° since there are no eclipses longer than 
about 0.03 in phase. (Light curve modeling is dis- 
cussed further below.) 

For an assumed value of AK, the measured val- 
ues of P, if wd and K^a, and our estimate of Mwd 
from log g, we compute Q and M to tai, and we infer 
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the orbital inclination i. This is shown in the lower 
panel of Figure | for M WD = 0.35 M , 0.40 Mq, and 
0.45 Mq. The upper panel of Figure || shows Mm 
vs. AK. If Mwd = 0.35 M Q , AK must be less than 
w 7.2 km s" 1 to allow i < 90°, and if i < 80° as sug- 
gested by the absence of eclipses, AK <J 4 km s _1 , 
much less than our estimate of AK w 21 km s _1 . For 
M WD = 0.4 M and i < 80°, allowed values of AK 
(< 14 km s _1 ) are still smaller than 21 km s _1 . In or- 
der to have AK ;> 21 km s -1 , M W d ^ 0.45 Mq is re- 
quired. Taken at face value, this "large" white dwarf 
mass favors the white dwarf models with a "thick" 
hydrogen envelope (Mh = 4 x 10~ 4 Mq, Benvenuto 
& Althaus 1998). 

For AK = 21 ± 2 km s" 1 and M W d = 0.45 ± 
0.05 Mq, the mass of the M star is M M = 0.28 ± 
0.05 M Q , i = 77 ± 7°, and Q = 0.62 ± 0.08, where we 
have propagated the errors on ifwD and Kn a . The 
mass of the M star is not too different from the mea- 
sured masses of the dM4e components of CM Dra- 
conis: M = 0.237 and 0.207 M , with errors of 4% 
(Lacy 1977). When M WD = 0.45 M Q , increasing AK 
by 5 km s" 1 reduces M M by 0.015 M . 

To summarize the argument: so far we have dis- 
cussed a mass estimate for the secondary star based 
on orbital dynamics. The two uncertain quantities 
in the analysis are Mwd, made uncertain by the un- 
known H envelope mass, and AK . Adopting a plau- 
sible value for AK (based partly on the radii of the 
components of CM Dra, which has a similar spectral 
type) leads to an estimate of Mm that is consistent 
with the more massive component of CM Dra. This 
is consistent, if slightly circular reasoning, but it sug- 
gests that the inclination of the binary orbit plane is 
fairly high. To avoid eclipses, however, i must not be 
too high, and thus higher stellar masses are favored, 
in order to account for the observed radial velocity 
amplitudes. 

An alternate method of estimating the mass uses 
an empirical mass-luminosity relation for M dwarfs, 
as follows. From the white dwarf's T c g and logg 
and assuming Mh = (Mwd = 0.40M Q ) we find 
log L wb /Lq = -0.5 (Althaus & Benvenuto 1997), 
and thus Mboi(WD) = 6.00. Using computed bolo- 
metric corrections and colors for hydrogen atmo- 
sphere white dwarfs from Bergeron et al. (1995; for 
T cff = 30,000 K and logg = 8), we then find 
My (WD) = 8.966 and M fl (WD) = 9.095. From the 
spectral decomposition described above, the M star 
contributes « 16% of the light in R, so Mr(M) = 



10.90. According to Kirkpatrick & McCarthy (1994), 
an M4.5V star has V — R = 1.37, so My (M) = 12.27. 
Finally, Henry & McCarthy (1993) give the relation 
logM M /M = -0.1681My + 1.4217 with an rms 
scatter in log M/M of 0.081. From this, M M = 
0.23 ± 0.04 Mq. The M star's mass inferred by this 
method is slightly higher (by sa 0.02Mq) if the white 
dwarf's hydrogen layer mass is Mh = 4 x 10~ 4 Mq, 
because the higher white dwarf mass (0.45M Q , Ben- 
venuto & Althaus 1998) implies a slightly larger ra- 
dius for the white dwarf (at fixed gravity), hence 
higher luminosities for both stars. Either value of 
Mm is consistent with the orbit-based mass estimate. 

As an aside, we can compute the distance to the 
system. The R magnitude of the system at phase 
0.0 (when the contribution of the heated face of the 
M star is minimal) is R = 16.22. The white dwarf 
contributes 84% of the light in R, hence the ap- 
parent R magnitude of the white dwarf by itself is 
m K (WD) = 16.22 - 2.5 log 0.84 = 16.41. The dis- 
tance is then d = 290 pc (0.4M Q white dwarf) or 
d = 308 pc (0.45M Q white dwarf). 

4.2. Light Curve Models 

The light curves described above show no obvious 
eclipses and are nearly sinusoidal. They are thus con- 
sistent with a very simple picture of the light modu- 
lation in PG 1224+309 being due to phase-dependent 
visibility of the bright, inward- facing hemisphere of 
a nearly spherical secondary star, which is uniformly 
emitting as a consequence of strong heating by the 
white dwarf. If one were sure of the radiative proper- 
ties of both illuminated and unilluminated sides of the 
M dwarf (i.e., if one knew the surface fluxes and limb- 
darkening coefficients in each color), one could model 
the light curves to infer constraints on the radius of 
the M dwarf and on the inclination. Such radiative 
quantities are becoming available for normal, unil- 
luminated photospheres of low-luminosity stars (cf. 
Claret 1998 for limb-darkening coefficients based on 
the Next-Gen model atmospheres of Hauschildt et al. 
1998). For the strongly illuminated hemisphere, how- 
ever, neither surface fluxes nor limb-darkening be- 
havior is known. Treating the problem by assum- 
ing that the atmosphere behaves "normally," i.e. with 
standard albedo and limb-darkening coefficients at an 
enhanced effective temperature, is not adequate for 
making robust conclusions regarding the secondary 
star (cf. Hilditch et al. 1996, who review problems of 
light-curve synthesis for strongly irradiated stars, in- 
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eluding the formal finding of limb- brightening in their 
analysis of the systems AA Dor and KV Vel). 

In view of this difficulty, we are content here to 
check whether the light curve amplitudes can be ap- 
proximately reproduced, using properties of the stars 
and the binary orbit consistent with our prior discus- 
sion, and without requiring such an extreme inclina- 
tion as to introduce an eclipse or other large distor- 
tion. Essentially, the question reduces to whether the 
pair of quantities i and Rm/o, that is needed to model 
the light curves is consistent with the properties of PG 
1224+309 as known from the discussion given above. 
For this purpose, a simplified treatment using black- 
body emissivities (at the local effective temperature) 
and a limb-darkening law appropriate for (unillumi- 
nated) cool photospheres may be adequate. 

We used the Wilson-Devinney (W-D; 1971) light 
curve synthesis code, fixing the temperature of the 
white dwarf at 29,300 K. The M star's temperature 
was fixed at 3100 K (unilluminatcd). The gravity 
darkening coefficients for the two stars were 1.0 and 
0.3, respectively. The bolometric albedos for reflec- 
tion heating and re-radiation were set to 1.0 and 0.5, 
respectively. For definitencss, a limb-darkening law 
was used that includes a logarithmic term and took 
the appropriate coefficients for the four filters from 
the tables in Van Hamme (1993). In the W-D code, 
the fractional radii R/a are expressed in terms of the 
mass ratio and the so-called ^-potentials of the two 
stars. The inclination and mass ratio were varied on 
a grid, 70° < i < 79° and 0.52 < Q < 0.71. For 
i > 79° eclipses are predicted. At each grid point, 
and assuming values for the fi-potentials, the B, V, 
R, and I light curves were fitted simultaneously, and 
the values of f2wD and Om were adjusted to find 
the minimum \ 2 - I n general, the best-fitting models 
underestimated the amplitude of the B light curve, 
matched more or less the amplitudes of the V and R 
light curves, and overestimated the amplitude of the 
/ light curve. In view of the modeling uncertainties 
in the case of strongly illuminated stars, the ability of 
the W-D code to roughly match the observed light- 
curve amplitudes for plausible values of i, Q and £Im 
is regarded as confirmation of the basic picture we 
have of PG 1224+309. 

4.3. Prospects for Improved Understanding 
of PG 1224+309 

Without the benefit of eclipses, how can accurate 
masses of the component stars of PG 1224+309 be 



derived? Additional observables include the M star's 
rotational velocity v mt sin i and a direct measurement 
of Km from absorption lines. The rotational velocity 
could be exploited, assuming synchronous rotation, to 
form the ratio Rm/o-m where cjm is the size of the M 
star's orbit around the center of mass. An absorption 
line velocity curve, depending on the lines used, would 
give an estimate of Km requiring a "K-correction" of 
the opposite sign to that used here. 

Additional information could be extracted from 
the light curves, given improved theoretical model- 
ing of irradiated atmospheres, to be incorporated into 
light-curve synthesis codes. In particular, believable 
light-curve amplitudes could be used to specify the 
radius of the M star, relative to the orbital separa- 
tion, as a function of the orbital inclination. Consis- 
tency between inferred mass and radius could then 
be demanded, to sharply limit the allowed range of 
inclination. 

In principle, one could measure the gravitational 
redshift of the white dwarf which would provide an 
independent constraint on the mass and radius of the 
white dwarf. The gravitational redshift would sim- 
ply be the difference in the two 7 values for the or- 
bital solutions of the two components. Currently, the 
measured 7 values are not known well enough, both 
in terms of formal statistical errors and in terms of 
systematic errors, to enable a meaningful measure- 
ment of the gravitational redshift (7wd = —4 ± 12 
km s _1 and 7H a = 10 + 6 km s^ 1 ; statistical errors 
are quoted). The Balmer lines have sharp emission 
features in their cores, which complicates the mea- 
surement of 7wd • Possibly the white dwarf has metal 
absorption lines in the ultraviolet, in which case a 
spaced-based measurement of 7wd might be possible. 

Finally, some improvement in the measurement of 
Kwd might be achieved from additional spectroscopy 
covering the entire orbital cycle. 

4.4. Evolutionary Considerations 

44.1. The Past: PG 1224+309 as a Post-Common 
Envelope Binary 

The present orbital separation of PG 1224+309 is 
a « 1.5 Rq. The present white dwarf was once the 
core of a giant star with a radius of « 100 Rq and 
therefore a > 100 Rq formerly. The standard picture 
of how such initially wide binaries can be drastically 
reduced in orbital size is the so-called "common enve- 
lope" (CE) process. Very roughly, in this picture the 
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envelope of the growing giant star engulfs the com- 
panion star, and orbital energy is lost due to friction 
against the envelope, leading to the possible ejection 
of the envelope. Iben & Livio (1993) present a well- 
organized and relatively recent review of the CE pro- 
cess and the stellar evolution scenarios and products 
involving it. Livio (1996) reviews how observations 
can further constrain the CE scenarios. 

The CE phase is very short-lived, and it unlikely 
to be observed directly, but post-CE binaries such as 
PG 1224+309 have the potential to illuminate some 
aspects of the CE process. One question of interest is 
whether a rule can be given for terminating the decay 
of the binary orbit — does the decay continue until 
essentially all material outside the (pre-white dwarf) 
"core" of the giant star has been expelled, and if so 
what is the mass of the residual hydrogen envelope? 
In PG 1224+309, M W d as inferred from T cff and log 5 
is ambiguous at the 0.05 Mq level, unless the hydro- 
gen envelope mass Mh is known. Given an estimate 
of Mwd, however, we can infer Mm along with the or- 
bital inclination, from AK and the mass ratio. With 
our best estimate of AK, the lack of eclipses favors 
a higher white dwarf mass, hence a relatively thick 
hydrogen envelope. PG 1224+309 misses being an 
eclipsing system by only a few degrees! Eclipses would 
provide additional constraints on the mass and radius 
of the white dwarf. It would then be possible to infer 
Mh with high confidence. This would be a new and 
valuable piece of information with which to confront 
common-envelope theories. 

Systems like PG 1224+309 can also sometimes be 
used to study the efficiency of the common envelope 
process. This is usually parameterized by the param- 
eter acE, which relates the initial and final orbits to 
the energy expended in ejecting the CE. The efficiency 
thus helps to determine the mapping of initial binary 
periods onto final (post-CE) periods. Except for the 
ambiguity in Mwd, PG1224+309 would be a useful 
system in potentially giving a lower limit to q;ce- The 
argument is as follows. The decay of the orbit must 
have released enough energy to overcome the binding 
energy of the envelope that formerly blanketed the 
white dwarf. The supply of orbital energy can be esti- 
mated from the present separation and masses (there 
is little dependence on the initial separation), and the 
present mass of the white dwarf gives the radius of its 
giant star progenitor at the moment that dynamical 
timescale mass transfer began (via the giant branch 
core-mass - luminosity relation). That radius, com- 



bined with an estimate of the envelope mass, gives 
the binding energy of the envelope. The ratio defines 
acE- £bind = a C EAE orhit or 

M 1 (M 1 -M 1R )_ ( M 1R M 2 M 1 M 2 \ 
Xa a r L aCE \ 2a f 2a a J 

using the definitions of Wcbbink (1984) as given in 
Livio (1996). Here Mi and M 2 are the initial masses, 
Mir is the core or remnant mass of the more evolved 
star, a and a/ are the original and final orbital sepa- 
rations, tl is the radius of the evolved star's original 
Roche lobe expressed as a fraction of a Q , and A = 0.5 
is a structure factor for that star's envelope. 

The calculated value of acE depends strongly on 
the assumed core mass Mm = Mwd- For PG 
1224+309, we first consider the case where Mm = 
Mwd = 0.40 M Q , M 2 = 0.25 Mq, and a f = 1.48%. 
Then the available orbital energy was about A£" or bit = 
3.4 x 10~ 2 in solar units with G = 1. Here we neglect 
the small second term on the R.H.S. of the displayed 
equation. For assumed initial masses of Mi = 1.0, 
1.2, and 1.4 Mq, the radius of the star when it con- 
tacted its Roche lobe and began mass transfer lead- 
ing to the common envelope was i?i bc = a a r^ = 75, 
71, or 68 Rq (using approximations from Tout et 
al. 1997). The binding energy of the envelope was 
approximately -Ebind = 1-6 x 10~ 2 , 2.7 x 10~ 2 , and 
4.1 x 10~ 2 , respectively. With these assumptions, the 
available orbital energy exceeded the binding energy 
only if Mi was less than 1.4 Mq, and even then most 
of the orbital energy was needed to unbind the en- 
velope; thus «ce cannot have been significantly less 
then unity. (Stars with values of Mi much smaller 
than those considered are less likely to have evolved 
so far yet. Also, for Mm = 0.4 Mq, which corre- 
sponds to a luminosity well below the tip of the red 
giant branch, mass loss on the giant branch is not 
likely to have been significant.) By this example we 
have shown how, in principle, observations of rela- 
tively wide post-CE binaries (with small negative or- 
bital energies) set a lower limit on «ce- However, in 
the specific case of PG 1224+309, the ambiguity in 
Mwd owing to our ignorance of Mh makes this par- 
ticular test less interesting. This is because increasing 
Mwd to w 0.45 Mq reduces the pre-CE binding en- 
ergy greatly and thus reduces the lower limit on «ce- 
Finding eclipsing, double-lined DA+dM systems may 
in favorable cases provide firmer constraints on acE- 
Analogous eclipsing systems would allow better con- 
straints to be placed on the component masses and 
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radii, which would in turn allow one to begin ad- 
dressing additional questions, such as whether post- 
CE secondaries satisfy a normal main sequence mass- 
radius relation. 

44.2. The Future: PG 1224+309 as a Pre- Cataclysmic 
Binary 

PG 1224+309 is a post-CE binary, as discussed 
above, and is also a pre-cataclysmic binary in the 
sense of Ritter (1986): the orbital separation a will 
shrink owing to gravitation wave radiation (GWR) 
until the M star contacts its Roche lobe and initi- 
ates mass transfer onto the white dwarf. (Angular 
momentum loss via magnetic braking, another pro- 
cess considered by Ritter, is thought not to oper- 
ate for low-mass M dwarfs.) Using formulae from 
Ritter (1986) and adopting M W d = 0.45 Mq and 
M M = 0.28 M , it will take about 9.0 x 10 9 years 
for GWR to bring PG 1224+309 into a semidetached 
configuration, at which time the orbital period will be 
2.5 hours. The present cooling age of the white dwarf 
is perhaps a few xlO 7 years (Althaus & Benvenuto 
1997), but when Rochc-lobc overflow (RLOF) mass 
transfer begins the white dwarf will have cooled to 
well below T cS = 5000 K. 

5. Summary 

We have measured the orbital period of PG 1224+309 
to be P = 0.258689 ± 0.000004 days. The semi- 
amplitude of the white dwarf's radial velocity curve is 
Kwd = 112+14 km s -1 . The velocity semi-amplitude 
of the Ha emission line arising from the irradiated 
hemisphere of the M star is Kn a = 160 ± 8 km s _1 . 
We estimate a "K-correction" of AK = 21 + 2 km 
s _1 , where Km = K^ a + AK. The implied mass 
ratio is then Q = M M /M WD = 0.62 + 0.08. The spec- 
tral type of the secondary is M4+. We fit synthetic 
spectra to the blue spectrum and find T e g — 29, 300 K 
and log g — 7.38 for the white dwarf. The white dwarf 
contributes 97% of the light at 4500 A and virtually 
all of the light blueward of 3800 A. 

The mass inferred for the white dwarf depends on 
the assumed mass of the thin residual hydrogen enve- 
lope: 0.40 < Mwd < 0.45 Mq for hydrogen envelope 
masses in the range < Mh < 4 x 10 _4 M q . Given 
the absence of eclipses, for likely values of AK and 
based on the observed values of i^wD and Km, we 
conclude that the mass of the white dwarf is closer to 
0.45 Mq. It thus appears that the white dwarf has a 



relatively large residual hydrogen envelope. The mass 
of the M star is M M = 0.28 ± 0.05 M Q if AK = 21 + 2 
km s -1 and M WD = 0.45 + 0.05 M & . In this case, the 
inclination is i = 77 ± 7°, only a few degrees away 
from giving rise to eclipses. Additional observational 
and theoretical work could improve our knowledge of 
the masses of the component stars. 

We argued that an accurately determined resid- 
ual hydrogen envelope mass can constrain theories 
of common envelope evolution, and we showed how 
accurate values of Mwd in binaries such as PG 
1224+309 can help set limits on «ce, the efficiency 
parameter used in common envelope scenarios. PG 
1224+309 itself will become a cataclysmic variable 
with an orbital period of about 2.5 hours, after a fur- 
ther « 10 10 years of orbital evolution. 
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